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What is claimed is : 

1. A method of imputing missing values in microarray 
data comprising the steps of: 

(a) clustering the data by a Gaussian mixture 
clustering model; and 

(b) estimating missing values by a GMCimpute 

i. ^orithm 

thereby imputing missing values in microarray data. 

2. The method of claim 1, wherein the Gaussian mixture 
clustering model comprises the steps of 

(a) determining a value of K; 

(b) partitioning the rows of the microarray data 
into K partitions; and 

(c) repeating a Classification Expectation- 
Maximization algorithm until the K partitions converge. 

3. A computer program product comprising a computer 
software program, wherein the computer software program, once 
executed by a computer processor, performs a method of 
imputing missing values in microarray data according to the 
method of claim 1. 

4. The computer program product of claim 3, wherein the 
Gaussian mixture clustering model comprises the steps of 

(a) determining a value of K; 

(b) partitioning the rows of the microarray data 
into K partitions; and 

(c) repeating a Classification Expectation- 
Maximization algorithm until the K partitions converge. 

5. A computer software program, wherein the computer 
software program, once executed by a computer processor, 
performs a method of imputing missing values in microarray 
data according to the method of claim 1. 
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6. The computer software program of claim 5, wherein 
the Gaussian mixture clustering model comprises the steps of 

(a) determining a value of K; 

(b) partitioning the rows of the microarray data 
into K partitions; and 

(c) repeating a Classification Expectation- 
Maximization algorithm until the K partitions converge. 

7. A computer comprising a computer memory having a 
computer software program stored therein, wherein the computer 
software program, once executed by a computer processor, 
performs a method of imputing missing values in microarray 
data according to the method of claim 1. 

8. The computer of claim 7 wherein the Gaussian mixture 
clustering model comprises steps of 

(a) determining a value of K; 

(b) partitioning the rows of the microarray data 
into K partitions; and 

(c) repeating a . Classification Expectation- 
Maximization algorithm until the partitions converge. 
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EM_ESTIMATE .... p. K , ^ S K> t, x K , A) 

{ FOR EACH ROW R OF /V WITH MISSING VALUES 

{ FOR/=1,...,K 

{ USE EM AND S,) TO ESTIMATE THE 

MISSING VALUES IN R 
Rj <- R WITH MISSING VALUES REPLACED BY ESTIMATES. 

} 

} 

R <- WEIGHTEDAVERAGE(f? 1 R K ). 

REPLACE R IN A BY R . 
RETURN A. 

} 

k_estimate(k; a) 

{ i* first part: initialization 7 

b *- rows of a without missing values. 

M"|> Mk> ^1 ^K> T 1' T K "* 

GAUSSIAN MIXTURE CLUSTERING OF B. 

A — EM.ESTIMATE (m ^ E, S K> Tl x K> A). 

I* SECOND PART: ITERATION 7 
REPEAT 

^.{.. \x. v .... 2^ S K , x v x K «— 

^ GAUSSIAN MIXTURE CLUSTERING OF A. 

' A EM.ESTIMATE fa, .... ^ L, S K t, x K , A'). 

} UNTIL CONVERGENCE 

} •: 

GMCimpute(S, A) 

{ FOR K= 1,2 S 

{ A K ^-K_ESTIMATE(K,v4). 
} 

RETURN (A,+ /4 2 + ... +A S )JS. 
} 
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